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Abstract 

Methods of Borel integration to calculate the binding ground energies and mass of bb quarkonia 
are presented. The methods take into account the leading infrared renormalon structure of the "soft" 
binding energy E(s) and of the quark pole masses m q , and the property that the contributions of these 
singularities in M(s) — 2m q + E(s) cancel. The resummation formalisms are applied to quantities which 
do not involve renormalon ambiguity, such as MS mass m q and a s (/x). 

1 Introduction 

The calculation of binding energies, masses of heavy quarkonia qq and another physical parameter using 
renormalon method has attracted the attention recently. The calculations, based on perturbative expansions, 
are primarily due to the knowledge of up to N 2 LO term (~af) of the static quark-antiquark potential V(r) 
^12] and partial knowledge of the N 3 LO term there, and the ultrasoft gluon contributions to a corresponding 
effective theory N 3 LO Hamiltonian PIEIH; and the knowledge of the pole mass m q up to order ~ a 3 @[2|. 
Another impetus in these calculations was given by the observation of the fact that the contributions of the 
leading infrared (IR) renormalon singularities (at b = 1/2) of the pole mass m q and of the static potential 
V(r) cancel in the sum 2m q + V(r) [S1IH1EI- Consequently, this cancellation effect must be present also 
in the total quarkonium mass M = 2m q + E qq ^2 or more precisely, in M(s) = 2m q + E(s) where 
E(s) is the hard+soft part of the binding energy, i.e., the part which includes the contributions of relative 
quark-antiquark momenta \k |, |k| > m q a s , i.e., soft/potential scales (predominant) and higher hard scales 
(smaller contributions). In addition, the binding energy has contribution E qq (us) from the ultrasoft momenta 
regime \k |, |k| ~ m q a 2 s . The ultrasoft contribution is not related to the b — 1/2 renormalon singularity, 
since this singularity has to do with the behavior of theory in the region which includes the hard (~ m q ) 
and soft/potential (~ m q a s ) scales. 

This contribution is based on ref. 1 1 3| and we present numerical calculation of the binding ground 
energies E b i (separately the s and the its parts) and the mass (2mb + E b i) of the heavy bb system, by taking 
into account the leading IR renormalon structure of rrn, and E bb (s), and combining some features as: (a) 
the mass that we use in the perturbation expansions is a renormalon- free mass ^3 El [01 El DS| which we 
choose to be the MS mass m q = fn, q (fi = m q ); (b) Borel integrations |12j are used to perform resummations. 
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Before resummations we perform separation of the soft/potential (s) and ultrasoft (us) part of the binding 
energies, and apply the renormalon-based Borel resummation only to the s part. The renormalization scales 
used in the Borel resummations are (ih ~ m q (hard scale) for 2m q , and m q a a < fi s < m q for E qq (s). 
The term corresponding to E qq (us) is evaluated at fx us ~ m q ot 2 s whenever perturbatively possible. Further, 
the Borel resummations are performed in three different ways: (a) using a slightly extended version of the 
full bifocal expansion Refs. ^J^]; (b) using a new "cr-regularized" full bifocal expansion introduced in 
the ref. ( c ) using the form of the Borel transform where the leading IR renormalon structure is a 

common factor of the transform ^1 EH (we call it R- method). The Borel integrations for both m q and 
E qq -(s) are performed by the same prescription (generalized principal value PV [201 EH EH 12] ) so as to 
ensure the numerical cancellation of the renormalon contributions in the sum 2m q + E qq (s). Furthermore, 
we demonstrate numerically that in the latter sum the residues at the renormalons are really consistent with 
the renormalon cancellation when a reasonable factorization scale parameter for the s-us separation is used, 
while they become incosistent with the aforementioned cancellation when no such separation is used. The 
numerical results allow us to extract the mass mj, from the known T(15) mass of the bb system. 



2 Pole mass 

Here we calculate the pole mass m q in terms of the MS renormalon- free mass rn q = m q (fx — m q ) and 
of a s (fi, MS), using elements of the renormalization group approach of Ref. and the bifocal expansion 
method Refs. ^|^]. The ratio m q /m q has, divergent, perturbation expansion in MS scheme which is at 
present known to order (Ref. [B] for ~a^; jjj for ~a 3 . Let us define 

S = ^.-l = U(fj,)[l + a( fJ ,)r 1 (n)+a 2 (n)r 2 (f,) + 0(a 3 )] , (1) 

ri(n) = Ki + P L m (fj,) , (2) 

t-2 (m) = K 2 + (2Ki/3o + f3i)L m (ii)+PlL 2 m (ii) , (3) 

(4/3)«i = 6.248/?o - 3.739 , (4) 

(4/3)k 2 = 23.497/3^ + 6.248/3i + 1.019/3 - 29.94 , (5) 

where L m = ln(/x 2 /m q ) , while (3 = (11 - 2n//3)/4 and ft = (102 - 38n//3)/16 are the renormalization 
scheme independent coefficients with nf = ni being the number of light active flavors (quarks with masses 
lighter than m q ). The natural renormalization scale here is fi = fih ~ m q (hard scale). 

Let us consider S as a function of the running coupling a(Q) — a(q)/w, and the perturbation expansion 
for S reads: 

oc 

S(a(Q)) = s «« n ' ( 6 ) 

Tl = 

which has to be summed. Then, one defines the Borel transform as 

B(b,Q) = y2s n -. (7) 

* — ' TV. 

n=0 

This series has a finite radius of convergence in the 6-plane, and we introduce the Borel function S(a) (Borel 
Integral) corresponding to S(a) as 

/•OO 

S(a) = / dbB(b,Q)e~ b/a (8) 
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and in our case Bs(b) is known to order ^b 2 



B s {b-a) = - 



r-2 



0(b 3 ) 



(9) 



It has renormalon singularities at b = 1/2,3/2,2,...,— 1,-2,... |14l 1221 12~3"| . The behavior of B$ near the 
leading IR renormalon singularity b = 1/2 is determined by the resulting renormalon ambiguity of m q which 
has to have the dimensions of energy and should be renormalization scale and scheme indepenent - the 
only such QCD scale being const x Aqcd EH- The Stevenson scale A [25] can be obtained in terms of the 
strong coupling parameter a(/i; c 2 , C3, . . .) = a s ([i; c%, C3, . . .)/tt, where Cj — 0j/0o (j > 2) are the parameters 
characterizing the renormalization scheme, by solving the renormalization group equation (RGE) |25| 



da(fi) 
din fi 2 



-0 o a 2 (p)(l + cia(n) + c 2 a 2 {u) H ) 



1 



W0J0 ^ 

1 



1 



1 



Cl 



/3 a(/i) O 



x 2 (l + c\x + C2X 2 + 
. 1 + cio(//)' 

In 



•) X 2 (\ + C\X) 



A = jj, cxp — 



1 



20 o a(a) 



1 + cxa{n)\ v 



exp 



1 

"20o 



o(ju) 



dx 



cia(n) 
(c 2 + C32; + c A x 2 H ) 



(1 + cix)(l + ax + c 2 x 2 H ) 

where v = a/(20o) = 0i/(20q). Expansion of expression l(T2"|) in powers of a{a) then gives 

1 



where 



r 3 



A = [i exp 



(4-C2) 



20oa(v) 



1 + r ka k (fJ.) 



k=l 



20o 



^2 = — 52 M - c 2) 2 - 2/?o(c? - 2 Cl c 2 + c 3 )] 



;[(c? - c 2 ) 3 - 60 a (cj - ca)(c? - 2cic 2 + c 3 ) 



48/3 31 

+8/3g(c£ - 3c?c 2 + C 2 + 2cic 3 - c 4 )] . 
The singular part of the Borel transform Bg{b) around 6=1/2 must have the form 



Bs{b\y) = N m ir-_ 



(l-2b) l + v 



l + £c fc (l -2b) k 



fe=i 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



and B ( s an ) (b; li) is analytic on the disk \b\ < 1. The MS coefficients for n/ = 4: a = 1.5400, c 2 = 3.0476 
and C3 = 15.0660 are already known [20 E], but for C4 we have only estimates [2B1GII] obtained by Pade- 
related methods. We consider C4 = 40 ± 60. Thus, Cj can be obtained: c\ = —0.1054, c 2 = 0.2736 and 
5 3 = 0.01 ±0.17. 

The bilocal method j^T| consists of taking in the expansion (|15fl for the analytic part B^ n '^ a poly- 
nomial in powers of b, so that the expansion of B$ around b = agrees with expansion ij^J). For that, the 
residue parameter N m in Eq. (|15|) has to be determined. Using Refs. |31| : 

N m = ^-R s (b= 1/2) , (17) 

and according to l|15|) 

R s (b- i i) = (l-2b) l+l 'Bs{b- l i) . (18) 
Then, the N m is estimated (see Ref. |U), using Rs(b) TPS and Pade approximation [1/1]: 

N m (n f =4) = 0.555 ±0.020 , (19) 

The bilocal expansion l|15|) has then for the analytic part the polynomial 

4^ /4 m) 
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*r - I--" (-) 

^ n— x 7 



where, by convention, ro = Co = 1. Then, the bilocal formula, is 
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Bs(M (bUoc ° = N m ir-^- L^- 

V m 9 (1 - 2&) 1+,y 



l + ^H fc (l-26) fe 



fc=i 



^ kl (3 k 1 ' 



fc=0 



Applying the (generalized) principal value (PV) prescription for the Borel integration 

»oo±ie / 1 



S(a) = — Rc / db exp 

Po ./±ie 



A) a Of) 



Bs(6;a*), (23) 



we obtain the pole mass m q in terms of the mass m q . The numerical integration is performed, using the 
Cauchy theorem ( Refs. [T^]). 

In Figs, n (a) we present the resulting (PV) pole masses of the b, as function of the renormalization 
scale \i. The spurious /^-dependence is very weak. In addition, results of another method ("R" -method ) are 
presented in Figs. Q (a), with the /^-dependence stronger in the low-/i region (fx/m q < 1). The R- method 
(Refs. ^3E3) cons i s ts in the Borel integration of the function l|18l) 

S = T^c f " * exp (>) «°ML , (24) 
A) J±ie \ foafa)) (l-25)i+- v ; 

where for Rs(b) the corresponding (NNLO) TPS is used. When we take rn\, = 4.23 GeV and we vary the 
values of the residue parameter N m , the bilocal method gives, at fJ,/m q — 1, variation Snib = T3 MeV. When 
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the central values of N m are used, the variation of the obtained values of m q with fi, when fi/m q grows from 
1.0 to 1.5, is about 5 MeV for mb (for i?-method: 4 MeV). When a is varied, the variation is about =f2 and 
=pl MeV for mb- The uncertainty in a s can be taken as a s (m T ) = 0.3254 ± 0.0125 32 , corresponding to 
a s (Mz) = 0.1192 ± 0.0015. This uncertainty is by far the major source in the variation of the pole masses: 
(8m b ) as =±\H MeV for bilocal method (±\% MeV for E-method). 

The natural renormalization scale \i here is a hard scale (j,~m q , and will be denoted later in this work 
as n m in order to distinguish if from the "soft" renormalization scale fj, used in the analogous renormalon- 
based resummations of the (hard+)soft binding energy E qq (s) (m q > /x > m q a s ) in Sec. 01 The fact that 
the two renormalization scales are different does not affect the mechanism of the (b = 1/2) renormalon 
cancellation in the bilocal calculations of the meson mass (2m q +E qq (s)), because the renormalon ambiguity 
in each of the two terms is renormalization scale independent ~ A, as seen by Eqs. (|13ll -l|2l. On the other 
hand, if i?-type methods (|24H [cf. also Eq. I|17|)] are applied for the resummations of 2m q and E qq (s), the 
renormalon ambiguities are renormalization scale independent in the approximation of the one-loop RGE 
running, and the renormalon cancellation is true at this one-loop level. 



3 Separation of the soft and ultrasoft contributions 

The perturbation expansion of the (hard + soft + ultrasoft) binding energy E qq of the qq heavy quarkonium 
vector (S = 1) or scalar (S — 0) ground state (n = 1, 1 = 0) up to the N 3 LO 0(m q a b s ) was given in [331 T ne 
reference mass scale used was the pole mass m q . The ground state energy expansion has the form 



4 r 
E qq = -g TO g 7r2a2 (/ i ){ 1 + [ fc i,o + h,iLp(lJ>)] 



4 
— i 
9 

+a 2 {p) [fc 2 , + k 2 sL p (p) + fe 2 ,2£p(/-0] 

+a 3 (v) [*3,o + hsL p (p) + h, 2 L 2 p (p) + k 3<3 L 3 p (fx)] + 0(a 4 )} (25) 

where 

LM = lnf 4 M . . ) • (26) 

The coefficients kij of perturbation expansion (|25H of the quarkonium (n = 1; £ = 0; S = 1 or 0) are given 
below. [Ml 03 EH 03: E3| - For N c = 3: 

/97 11 \ 

fei,i = 4/3 , fei,o = ( y - Y n f ) ' ( 27 ) 

927 193 

k 2 ,2 = 12/3g , k 2A = — —n f + n}, (28) 

fc 2j0 = 361.342 - 40.9649 n f + 1.16286 n) - 11.6973 S(S + 1) , (29) 

„o , 4521 10955 1027 2 5 , 
fe 3 , 3 = 32/3 3 , %2 = ^ 24~n/ + —n) - -n) , (30) 

fc 3) i = 7242.3 - 1243.95 n f + 69.1066 - 1.21714 n) 

+ 25^2 (_67584 + 4096 n/) S[S + 1] ' (31) 
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(7839.82 - 1223.68 n f + 69.4508 n) - 1.21475 n)) 
-(-109.05 + 4.06858 n f ) S(S + 1) 

2 

~ (-1089 + 112 S(S + 1)) In «/*)) + 2^ 



(32) 



Here, a3 have been estimated in Ref. |31j . obtained from the condition of renormalon cancellation in the sum 
(2m q + V qq (r)) 

— a 3 (n f =4) «86. ± 23. . (33) 

The coefficients fcij in the expansion (|25|l originate from quantum effects from various scale regimes of the 
participating particles: (a) the hard scales (~m q ); (b) the soft and potential scales where the three momenta 
are |q| ~m q a s (\q°\ ^m q oi s in the soft and \q°\ ~m q a 2 s in the potential regime); (c) ultrasoft scales where \q°\ 
and |q| are both ~m 9 a 2 . The coefficients are dominated by the soft scales; the hard scales start contributing 
at the NNLO |3] and are numerically smaller. For this reason, we will usually refer to the combined soft and 
hard regime contributions to the binding energy as simply soft (s) contribution E qq (s). Strictly speaking, 
it is only the pure soft regime that contributes to the 6=1/2 renormalon. However, for simplicity, we 
will resum the hard+soft contributions E qq (s) together. This will pose no problem, since the hard regime, 
being clearly perturbative, is not expected to deteriorate the convergence properties of the series for E qq (s). 
The natural renormalization scale /x in the resummations of E qq (s) is expected to be closer to the soft scale 
(m q < a s (i < m q ). 

On the other hand, the N 3 LO coefficient k^^ obtains additional contributions from the from the 
ultrasoft (us) regime. The leading ultrasoft contribution comes from the exchange of an ultrasoft gluon in 
the heavy quarkonium (HUHj- It consists of two parts: 

1. The retarded part, which cannot be interpreted in terms of an instantaneous interaction 

1 . - , 2 /4 N 2 



-3fc 3)0 (^,ret.) = -— ( - I if w +41.014, (34) 

where if w -81.538 is the QCD Bethe logarithm - see Refs. 0E1. 

2. The non-retarded part can be calculated as expectation value of the us effective Hamiltonian TL US in 
the Coulomb (i.e., leading order) ground state |1), where 7i us (in momentum space) was derived in 
Refs. Direct calculation of the expectation value, here in coordinate space, then gives: 

1 9 1 2 
^■fc 3i o(ws,nonret.) = --rs; (l\H us \l) 



2]x{- 2 ^fa 3 M(l\ 1 -\l) 



47I" 5 m q a 5 (/J,) 7T 5 TOqa 4 (/x) 
X i 4 , 5 , „i f 277r3 3/ wn 1 

-17^<l|il> + *£M mr )\l) 
m q r z 3 m~ 

+37r^(l|{A r ,i}|l)} , (35) 



= -14.196 



In ( ^r—- ) +0.9511 

m q aj(u) 



(36) 
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Here, = — (4/9)m g a 2 (/i) is the Coulomb energy of the state |1), and /if is the factorization energy 
between the soft (~ m q a s ) and ultrasoft (~ m q a 2 s ) scale. 

The s-us factorization scale [if can be estimated as being roughly in the middle between the s and 
us energies on the logarithmic scale [2] 

= k m q a s (n s ) 3/2 , (37) 

where k ~ 1 and fi s « i?s Therefore, the ultrasoft part of the N 3 LO coefficient k 3y Q can be rewritten, 
by Eqs. (|34|) . J5J and i|37|) . in terms of the s-us parameter k as 

-^k 3Q (us) = 27.512 + 7.098 ln(a s (fj, s )) - 14.196 ln(/t) . (38) 

The soft scale /z s appearing here will be fixed by the condition /i s = (A/3)m q a s (p s ). 

The formal perturbation expansions for the separate soft and ultrasoft parts of the ground state 
binding energy 1)25(1 are then 

4 2 * 

E qq {s) = --m g Tr 2 a 2 (iJ,){l + J2 ai (l J -)J2 k ^ L P^y 

i=l 3 =0 

3 

+a 3 (n) fc 3,, + a 3 (M) [*3,o-*3,o(u*)] + 0(a 4 )} , (39) 

3=1 

£ g? (us) = --m g 7r 2 a 2 ( A i){a 3 (Ai)fc 3 ,o(ws) + e'(a 4 )} . (40) 
9 

The energy E qq (s) (|39|l contains the leading IR renormalon effects, and E qq {us) (|40|l does not. In these 
expressions, the common factor is the soft scale fi p (p) = (4/3)m q a s (/z) which is also present as the reference 
scale in the logarithms L p (n) = ln(/z//Lt p (/x)) appearing with the coefficients kij (when j > 1) in Eqs. i|25|) . 

We will re-express m g everywhere in with the renormalon- free mass fn q , and will consider the 
dimensionless soft-energy quantity E qq (s)/m q . 

Thus, we will divide the soft binding energy with the quantity ~p(ji) = (4/ '3)m q a s (jl) , where fi can 
be any soft scale. We will fix this scale by the condition Jl — (A/3)rn q a s (]l) (^> = fi s ). Further, in the 
logarithms L p (/j,) we express the pole mass m q in terms of m q and powers of a(/i) (cf. SecEJ), and the powers 
of logarithms ln fc [a(^)] we re-express in terms of ln fe [a(/I)]. This then results in the following soft binding 
energy quantity F (s) to be resummed 

F(s) ee -^Jki^l = °M t 1 + a 0*)/i + « 2 (m)/2 + a 3 {fi)f 3 + C(a 4 )] , (41) 
47T Z m q a{n) 

where the coefficients fj depend on In a(ji) and on three scales: the renormalization scale \i (> m q a s ), the 
(fixed) soft scale jtz, and m q . The coefficient f 3 depends, in addition, on the parameters k l(37|l - l(38|l . /x s , and 
03 ©■ The coefficients fj are written explicitly in Appendix 1X1 The 6=1/2 renormalon in the quantity 
F(s) is then of the type of the renormalon of the pole mass m q discussed in the previous Sec. |5] 

However, if we divided in Eq. I|41|) by m q instead of rn q and at the same time used in the resulting 
fj -coefficients lnm g , the numerical resummations of F(s) by methods of Sec. 0] would give us values for 



/'./■ 



(Eg Ei 



us 



,1/2 
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E q q(s) different usually by not more than O(10 1 MeV) (we checked this numerically). We will briefly refer 
to these approaches later in this Section as "pole mass" approaches. A version of such pole mass bilocal 
approach was applied in Ref. |12| for resummation of the unseparated E qq (s+us). 

The ultrasoft part (|40|) . on the other hand, has no b = 1/2 renormalon. The mass scale used there 
should also be renormalon free (nig). The renormalization scale fi there should be adjusted downward to the 
typical us scale of the associated process \i >— > fi us (^m q a 2 s ) in order to come closer to a realistic estimate 

E qq -(us) w -^m 9 7r 2 /c 3i0 (us)a 5 (/i tis ) . (42) 



4 Mass rrib from the known mass of the vector bb [T(1S)] 

The soft binding energy quantity to be resummed is F(s) of Eq. (|41|l . However, in the N 3 LO coefficient /3 
we have dependence on 03, and on the s-us factorization scale parameter jc~l. Then, the value of k (|37ll 
is obtained requiring that the residue parameter values be reproduced from the Borel transform of the soft 
binding energy quantity F(s) of Eq. l|4"T|) . 
Similarly as in Eq. I)15|l . we have 

+ 4? S )(M («) 

where the factor in front of the singular part was determined by the condition of renormalon cancellation of 
the sum 2m q + E qq (s). We now define in analogy with Eq. I|18|) 

R Fis) (b;n;fx f ) = (l-2b) 1+l/ B F{s) (b-^ t i f ) . (44) 

Here we denoted, explicitly the dependence on the factorization scale ///■ and 

2tt m n a(71) „ . , , , 

m = T^T^W^w)!^ • (45) 

Theoretically, Rp^ s )(b) should be a function with only a weak singularity (cut) at b — 1/2, and the nearest 
pole at b — 3/2 ( |2E|)- Resummations such as Pade approximations (PA's) are applied, then [2/1] (6) 

has physically acceptable pole structure |6 po ic| > 1 for most of the values of /1 > m q a s and K ~ 1. and the 
residue parameter N m is reasonably stable under the variation of /x. 

In Fig. Eta) we show the dependence of N m on k, at a typical /j, value /x = 3 GeV, for the bb system. 
The know value of N m is obtained by the Rf( s ) [2/l](6 = 1/2) expression at k w 0.59. In Fig. Elb) we 
present, for k = 0.59, the dependence of calculated N m on the renormalization scale /i. There, we include 
also the ([2/l]-resummed) curve for the case when no separation of the s and us parts of the energy is 
performed. In that case, the obtained values of N m are unacceptable. The other values of the input 
parameters are chosen to have the bb "central" values: a 3 /4 3 = 86; mj = 4.23 GeV; Jl = 1.825 GeV (« tt s ) 
and a s (Jl;nf = A) = 0.3263(« a s (/x s ;n/ = 4) = 0.326) [from: a s (m T ; n/ = 3) = 0.3254, i.e., a s (M z ) = 0.1192 
1221 . For the RGE running, we use four-loop MS /3-function (TPS). 

Variation N m — 0.555 ± 0.020 [bb] implies k — 0.59 ± 0.19. If, on the other hand, 03 parameter is 
varied, then for bb k = 0.59 =F 0.06. Thus, the value of s-us factorization scale parameter n is influenced 
largely by the allowed values of the renormalon residue parameter, and significantly less by the allowed values 



2tt m q a(jl) (l-2b) 1+ " 



l + ^? fc (l-26) fc 



fc=i 
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a3 of the N 3 LO coefficient of the static qq potential. Therefore, we will consider the variations of N m and of 
k to be related by a one-to-one relation, while the variations of 03 will be considered as independent. 
In this way, we have the following value for the s-us factorization scale parameter: 

N m = 0.555 ± 0.020 k = 0.59 ±0.19 (ri/ = 4, 5 = 1) (46) 

and thus we obtain the N 3 LO TPS l|41|l for the soft part of the ground binding energy. Now with the value 
of k, we can perform the resummation of the soft part of the ground binding energy. The full bilocal method 
[THEI can be performed as in Sec.|3 Eqs. and Therefore 

^r' <m - »-h *£M (!-»>». |> (1 - 2t)t + 1 fit » ■ <«' 

/c — A, — 

where the coefficients hk in the expansion of the analytic part are known up to order k = 3 

3 

hk = fk N m f (2/3 ) fc ^ gn r w + !l 1 =e (fc = 0,l,2,3). (48) 

The result have spurious /i-dependence, and for the of the Borel transformation turns out to have a problem- 
atic behavior: the obtained pole of B^? ') [2/1] (b) is unacceptably small in size: |6 po ie| < 1/2. Theoretically, 

B^ g )(b) should have the nearest pole at b — 3/2 |2B]- The reason for this problem appears to lie in the 
specific truncated form of the singular part taken in the bilocal method. While the latter part describes well 
the behavior of the transform near b = 1/2, it influences apparently strongly the coefficients hk and thus the 
analytic part, so that no reliable resummation of that part (apart from TPS) can be done. This problem 
can be alleviated by introducing a "form" factor which suppresses that part away from b * 1/2, but keeps it 
unchanged at b ~ 1/2. We choose a Gaussian type of function, and the "cr-regularized" bilocal expressions 
for the Borel transform is 



x exp 



£ wa " rV ■ ,49) 



fc=0 



The corrective terms 1 / (8a 2 ) and Ci/(8cr 2 ) in the coefficients of Eq. (|49|l appear to ensure the correct known 
behavior of the renormalon part up to order ~ (1 — 2b)~ u+2 . Numerical analysis indicate that a is between 
zero and one. Namely, when a decreases from <r — oo to about a w 0.3-0.4, the value of the pole of the [2/1] 
Pade-resummed analytic part Bf?f\b) of Eq. l(4^|) gradually turns acceptable (|6 po ie| > 1) an d rather stable 
when the renormalization scale /i varies in the interval [m q a s ,mq] (except close to \x ss m q a s ). When the 
value of a falls below 0.3, the analytic part starts showing erratic behavior again and the Borel resummation 
gives significantly differing results with the TPS- and the Pade-evaluated analytic parts. On these grounds, 
the obtained optimal a turn out to be 

cr = 0.36 ±0.03 (ti/ = 4, 5 = 1). (50) 
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In Fig.[3ta) we present the Borel-resummed soft part of ground state energy for the bottonium (S = 1), 
as a function of the a parameter of method (|49ll . 

Finally, the results for the soft binding energy E bb (s) of the ground state of bottonium using the R- 
method |18U19) . where we resum the function R F ^(b; /i) (|44|l and then employ the (PV) Borel resummation 
as written in Eq. (|24[1 (with i£pv s ) instead of Rs there)., as functions of the renormalization scale u, are 
presented in Fig.QIa). We observe from the Figure that the bilocal "cr-regularized" method (|49|l (a = 0.36) 
gives the TPS and PA results closer to each other. The methods er-TPS, er-PA, and i?-PA give similar results 
in the entire presented /i-interval. i?-TPS appears to fail at low /i (« m\,a s ~ 1-2 GeV). In Fig. 0^b) we 
include, for comparison, the simple TPS evaluation of E bb (s), according to formula [cf. Eq. 141|) ] 

F(s) (TPS) s _ 9 = J__ {s) = a{ ) [x + a(fi)fi + a 2 [ )h + o 3 ( )/a ] (51) 

47r m b a s (fi) 

where for N 2 LO TPS case we take f% = 0. In Fig. 2Jb) the same input parameters are used as in Fig.^Ja). 
We see that the perturbation series shows strongly divergent behavior already at N 3 LO. In this Figure, we 
also included the "perturbative" ultrasoft part E^\us;fj,) calculated according to [see Eqs. (|38|) and 140|) ]. 

F (p) M = _ 9 _J_ B (p) (ua;/i ) = fc 3oa 4 (/i) . (52 ) 
4.1 Extraction of bottom mass 

The estimate of the perturbative part is given in Eq. I|42|l . where it was essential to take for the renormal- 
ization scale a us scale \i ~ fi us ~ m q a 2 s . 

For the bottonium case, this scale is below 1 GeV, the energy at which we cannot determine pertur- 
batively a s (p,). This indicates that in the bottonium the us part of the binding energy has an appreciable 
nonperturbative part. The lowest energy at which we can still determine perturbatively a s is /i « 1.5-2.0 
GeV, giving a s (/z) ~ 0.30 — 0.35. Although this is a soft scale for bb, we will use this also as an ultrasoft 
scale. Then by Eq. |@2) 

E b - b (us)^ « -^m q TT 2 k 3 , (us)a 5 (nus) ~ (-150 ± 100) MeV . (53) 
The nonperturbative contribution coming from the gluonic condensate is given by 

E^us)^ « m b 7r 2 ^ (^m b a s ( Pus ^j (a^G^G^) » (50 ± 35) MeV , (54) 

where we used m& = 4.2 GeV, and the value of the gluon condensate ((a s /n)G 2 ) = 0.009 ±0.007 GeV 4 |4T1] . 
Eqs. |J5j5J and (JSJl give 

E bb {us) (p+np) w (-100 ± 106) MeV , (55) 

The finite charm mass contributions has been calculated in Ref. (Refs. |Sl l42ll4^] '). The contribution to 
the mass Mr (IS) — (2m b + E bb ) is 

5M T (1S, m c 0) w 25 ± 10 MeV . (56) 

The estimates (|55|l , and (|56|l then give a rough estimate of the us and m c ^ contributions to the bottonium 
mass 6M~c(lS; us + m c ) « (—75 ± 106) MeV. The mass of the T(IS') vector bottonium ground state is well 
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measured My(lS) = 9460 MeV with virtually no uncertainty 03]. Therefore, the pure perturbative "soft" 
mass is 

M T (IS; s) = 2m b + E bb (s) = 9535 T 106 MeV , (57) 

where the uncertainty is dominated by the uncertainty of the us regime contribution. Our numerical results 
for E bb (s) and for m& allow us, by varying the input value of mj, to adjust the sum 2m b + E bb (s) to the 
value given in Eq. 1|57[) . For the soft binding energy we apply the "cr- regularized bilocal methods er-TPS and 
cr-PA, and R-TPS and R-PA, with the input parameters: a s (M z ) = 0.1192; ju = 1.825 GeV_(« fi s ), thus 
a s (/I,n/=4) = 0.3263 [a s (fj, s ) = 0.326]; N m = 0.555; k = 0.59; a = 0.36; a 3 /4 3 = 86; c 4 (MS) = 40. For 
2m b we apply the bilocal-TPS and i?-TPS method, with renormalization scale /i m /fS(, = 1, both methods 
giving us very similar results [cf. Fig. [I] (a)]. The bilocal-TPS method is applied for 2m b when cr-TPS and 
cr-PA are applied for E bb (s); the -R-TPS is applied for 2m b when .R-TPS and R-PA are applied for E bb (s). 
The extracted values of m b = m b (fj, = rn b ) are then 



m b = 4.225 ± 0.054 GeV 


(cr-TPS) , 


(58) 


= 4.220 ± 0.056 GeV 


(Cr-PA) , 


(59) 


= 4.243 ± 0.080 GeV 


(.R-TPS) , 


(60) 


= 4.235 ± 0.068 GeV 


{R-PA) . 


(61) 



Finally, the average of the central values gives us 

m b = 4.231 ± 0.068 GeV . (62) 

5 Comparisons and conclusions 

Our result for the mass m b = 4.231 ± 0.068 GeV, is compared in Table with recently values of m b . 
The only input parameter common to all these methods is a s . Our central value was a s (m T ) = 0.3254 
a s (Mz) — 0.1192] since such [35], or similar 05] ^5], values follow from the (nonstrange) semihadronic 
t decay data which are very precise 0B]. On the other hand, the world average as of September 2002 is 
a s (Mz) = 0.1183 ± 0.0027 [17] . Most of the authors during the last four years used central value a s (Mz) ~ 
0.118. Therefore, for comparisons, we convert our results H58JI to this central value of a s - more specifically, 
from a s (M z ) = 0.1192 ± 0.0015 to 0.1180 ± 0.0015. This can be easily done by inspecting in Table[T]thc 
column under a s , an increase in the central values of m b by 11, 12, 8 and 10 MeV, respectively. This gives 
the average 10 MeV higher than in Eq. (|62|l 

m b = 4.241 ± 0.068 GeV average when: a s {M z ) = 0.1180 ± 0.0015 . (63) 

There are two important numerical effects in our result. The first is the separate evaluation of the 
"perturbative" ultrasoft energy part at the corresponding low renormalization energy (1.5 — 2 GeV), Eqs. (|42|l 
and l|53(l . If we had not separated the ( "perturbative" ) ultrasoft from the soft part of the binding energy, 
the use of the common renormalization energy scale fj, (« 3 GeV) in the resummation then would have given 
us the central value of E bb (us) by about +100 MeV higher. Then the extracted value of rn b would have gone 
down by about 46 MeV, giving the value rn b « 4.195 ± 0.068, with the central value close to that of L03 in 
Tabled On the other hand, that latter value is quite clearly lower than the value PS02 in Tabled by about 
150 MeV, principally because of the 6=1/2 renormalon effect which were taken into account here and in 
Ref. ^21- Thus, the renormalon effect brings down the extracted central value of m b by about 150 MeV, but 
the separate evaluation of the ultrasoft contribution brings it up by about 50 MeV. 
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A Coefficients for the expansion of the soft binding energy 

We write down here the explicit coefficients fj of the expansion (|41|l for the soft part of the ground state 
binding energy. The logarithms appearing in these expressions involve three scales [/i, Jl, rn q and ~jl(Ji) = 

(4:/3)m q a s (Ji)] 



In = ln[=^J , X 2 = m^f) , L„ = ]n[^) ■ (64) 



The coefficients fj are 



fi = i(35 + 22L 1 -HL M -llL 2 ) + i(-ll-6Li + 3L Al + 3L 2 )n / . (65) 



h = f^+f^nf+f^nj, (66) 
/ 2 (0) = [381.67 + 90.75L? + 30.25L^ + l x (246.42 - 121L M - 60.5L 2 ) - 48.5L 2 

+L M (-205.25 + 60.5L 2 ) - 11. 6975(5 + 1)] , 
ffi = [- 42.7469 -11L?-3.6667L^ + L M (26.6944-7.3333L 2 )+6.8056L 2 

+Li(-33.0556 + 14.6667L P + 7.3333L 2 )] , 
fi 2) = [1.16286 +(3/9)l;+(1/9)LJ + Li(1-(4/9)L p -(2/9)L 2 ) 

+£ M (-0.81482 + (2/9)L 2 ) - 0.18518L 2 ] . 



h = f^ + fPnf + ftnj + f^nj, (67) 
/ 3 (0) = [6726.11 + 665.5^- 166.3757^ (40.802+ (-10.599 + L^L^) 
+L?(2381.5 - 1497.38^ - 499.125L 2 ) - 871.429L 2 
— 499. 125(— 1.884+ L^)L ll L 2 - 201.438^ + L 1 (7457.17 
-497.292i 2 + £ M (-4346.38 + 998.25L M + 998.25L 2 )) 
-257.341(0.2112 + L X - 0.75L,, - 0.25L 2 )S(S + 1) 
-61.4109 (-6.1394 + S{S + 1)) In {a s {fj, s )) + 440.172 ln(/t) + 2a 3 /4 3 ] , 

/ 3 (1) = [- 1274.33- 1277.92Li -471. Ylhh\ - Yl\h\ + 1182.32^ 
+843.667iiL M + 272.25L^ t - 335.813L^ - 181.5LiL^ 
+30.25X 3 + 124.501L 2 + 108.361LiL 2 + 90.75if£ 2 
-186.708X^2 - 181.5X1^2 + 90.75£^L 2 + 36.729^ 
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+ (4.0686 + 15.5964Li - 11.6973^ - 3.8991L 2 )S(S + 1)] , 
/| 2) = [70.8992 + 70.2453Li + 28.9722L? + 7.3333L? - 65.9925L M 

-51.6667LiL M - 16.5L?L M + 20.5972^ + WL-J?^ _ 1.83331^ 
-5.1939L 2 ~ 6.5741Lii 2 - h.hL\h 2 + 10.9167L M L 2 
+ll J LiL AI L 2 - h.hL\L 2 - 2.0972L|] , 

/ 3 (3) = [- 1.21475 - 1.21714L X - (5/9)L? - 0.14815L? + 1.16286^ 
-\-L\L^ + (l/3)i?i M - 0.407412,2 _ (2/9)£ii£ + 0.03704LJ 
+0.05429L 2 + (l/9)iii 2 + (l/9)L?L a 
-0.18518X^2 - (2/9)ii£ M i 2 + (l/9)£*£ a + 0.03704^] . 
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biloc.(N m =0.535) 
biloc.(N m =0.555) 
biloc.(N m =0.575) 
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> 

a> 

O 4.95 



m b = 4.23 GeV 
a s (M z )=0.1192 



4.94 



4.93 



4.92 



0.6 
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Fi gurG 1: The (PV) pole mass of the bottom quark, as function of the renormalization scale [i. The input parameters used were 
rfTj = 4.23 GeV and the residue parameter values Eq. 1191 . The reference value for a s (in MS) was taken to be o 3 (m T ) = 0.3254 
( 1321 1 corresponding to a s (M z ) = 0.1192. 
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Figure 2: The residue parameter value N m as calculated from the soft part of the binding energy of the 
bottonium according to Eq. (|45|l . (a) as a function of the s-us factorization scale parameter n (|37[l . at fi = 3 
GeV; (b) as a function of the renormalization scale fi, at k = 0.59. Further explanations given in the text. 
In Fig. (a), the known values of N m are denoted as dotted horizontal lines. 
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Figure 3: (a) Soft part of the ground state binding energy of 66, evaluated with the (PV) Borel-resummed 
expression (|49|l . as a function of the method parameter a. Details are given in the text. 
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Figure 4: (a) Soft part E b i(s) of the ground state binding energy of bb, evaluated with four different methods 
involving (PV) Borel resummation, as functions of the renormalization scale /i. Details are given in the text. 
In Fig. (b) the simple TPS results for E bb (s) are included [Eq. l|51[l ]. as well as the "perturbative" ultrasoft 
part E^f (us; (j.) [Eq. 



Table 1: The separate uncertainties 8m b (in MeV) for the extracted value of uib from various sources: 1.) us 
[6E bb (us)^ +n ^ = -100±106MeV]; 2.) ^ = 3±lGeV; 3.) fi m = m b (l ±0.5); 4.) a a (m T ) = 0.3254 ±0.0125 
[a,(M z ) = 0.1192±0.0015]; 5.) N m = 0.555±0.020 [k = 0.59±0.19]; 6.) a 3 /4 3 = 86. ±23.; 7.) c 4 = 40.±60.; 
8. cr = 0.36 ± 0.03; 9. m c (6M T (m c ^0) = ±10 MeV). 
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±45 
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Table 2: Recently obtained values of (MS) m& mass obtained from T sum rules or from spectrum of the 
T(1S) resonance. Wherever needed f |33l I12| 1. the central mass values were adjusted to the common input 
central value a s (Mz) = 0.118. 



reference 


method 


order 


m b (m b ) (GeV) 


PP98 EH 


T sum rules 


NNLO 


4.21 ±0.11 


MY98 01] 


T sum rules 


NNLO 


4.20 ±0.10 


BS99 HBJ 


T sum rules 


NNLO 


4.25 ±0.08 


H00 gSj 


T sum rules 


NNLO 


4.17 ±0.05 


KS01 [HI 


T sum rules 


NNLO 


4.209 ± 0.050 


CH02 E2]| 


T sum rules 


NNLO 


4.20 ±0.09 


E02 [FT] 


T sum rules 


NNLO 


4.24 ±0.10 


P01 mi 


spectrum, T(1S) 


NNLO 


4.210 ±0.090 ±0.025 


BSV01 |H] 


spectrum, T(1S) 


NNLO 


4.190 ±0.020 ±0.025 


PS02 [Sj 


spectrum, T(1S) 


N 3 LO 


4.349 ±0.070 


L03 [E] 


spectrum, T(1S) 


N 3 LO 


4.19 ±0.04 


this work, Eq. 


spectrum, T(1S) 


N 3 LO 


4.241 ±0.070 
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